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Abstract 

Define the nuclear dependence coefficient a{A, qx) in terms of ratio of trans- 
verse momentum spectrum in hadron-nucleus and in hadron-nucleon colli- 
sions: rf % / da . 2 = A a ( A ' qT > . We argue that in small qx region, the a(A,qx) 

for the Drell-Yan and J/i/j production is given by a universal function: a+bq T , 
where parameters a and b are completely determined by either calculable 
quantities or independently measurable physical observables. We demonstrate 
that this universal function a(A, qx) is insensitive to the A for normal nu- 
clear targets. For a color deconfined nuclear medium, the a(A, qx) becomes 
strongly dependent on the A. We also show that our a(A, qx) for the Drell- 
Yan process is naturally linked to perturbatively calculated a(A, qx) at large 
qx without any free parameters, and the a(A, qx) is consistent with E772 data 
for all qx- 
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I. INTRODUCTION 



It has long been observed that in high energy hadron-nucleus collisions, the transverse 
momentum spectra of produced particles differ significantly from that in hadron-nucleon 
collisions [l],@]. This anomalous nuclear dependence is known as the Cronin effect. In recent 
years, many data on such nuclear dependence for the Drell-Yan |3|-|5| and J/ip produc- 
tion [|]-§| became available, and these new data have renewed our interests to understand 
the observed novel effect |9|]. Although it is believed that the Cronin effect is a result of 
parton multiple scattering inside nuclear medium, systematic calculations in QCD for such 
anomalous nuclear dependence only exist for the large transverse momentum (q-r) region 
IT)|-|T2|, but not for the small qx region due to technical difficulties in handling multi-scale 



calculations in QCD perturbation theory. In this paper, we derive the nuclear dependence 
in transverse momentum distributions for the Drell-Yan and J/ip production in the small 
qx region by combining constraints from perturbatively calculable quantities in the large q? 



region and the available information from independently measured physical observables [13 



In high energy nuclear collisions, in addition to short- distance single scattering, parton 



multiple scattering becomes very important. According to QCD factorization theorem [|T4 
physical observables due to single hard scattering in hadronic collisions depend on the par- 
ton distributions. In terms of generalized QCD factorization theorem contributions 



to physical observables from parton multiple scattering are directly proportional to multi- 
parton correlation functions. These correlation functions are as fundamental as the parton 
distributions, and they provide complementary information on nonperturbative QCD dy- 
namics fl6|]. Since single hard scattering is localized in position space, it is the multiple 
scattering that is sensitive to the size of nuclear medium, and is therefore responsible for 
the anomalous nuclear dependence. Measurements of such anomalous nuclear dependence 
provide a window of opportunities to explore the dynamics of parton correlations, which 
have important implications for the physics in relativistic heavy ion collisions. 

The nuclear dependence of the transverse momentum spectrum is often presented in 
terms of the ratio of cross sections R or the nuclear dependence coefficient a. For example, 
for the Drell-Yan pair production, the R and a are defined by 

where Q and qx are the Drell-Yan pair's total invariant mass and transverse momentum, 
respectively, and A is the atomic weight of the nuclear target. In Eq. ([T]), da hA /dQ 2 dq T and 
da hN j dQ 2 dq T are the transverse momentum spectrum in hadron-nucleus and hadron-nucleon 
collisions, respectively. Since single hard scattering is localized in space, both R(A, qx) and 
a(A, qx) should be very close to one if there is no multiple scattering. However, data on 
the ratio of the Drell-Yan transverse momentum spectrum show that R(A, qx) very much 
differs from one (or a (A, q?) is significant away from one) In Ref. ||], data on R(A, q?) 
shows a nontrivial dependence on both A and q? for different nuclear targets (including C, 
Ca, Fe, and W). In terms of a (A, q?), the Drell-Yan data has the following general features: 
a(A, q T = 0) is less than one, it increases as qT increases, and it can be as large as 1.07 

When q T is large (~ Q), the g-r-dependence of a(A,q T ) can be calculated within QCD 



perturbation theory, and a significant nuclear enhancement in a(A, qr) was predicted [11 . 
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Although it was argued in Ref. [II] that in the small qx region, the a(A, qx) or the transverse 
momentum spectrum should show nuclear suppression, no quantitative analysis or prediction 
was given for the suppression. Since almost all existing data on the nuclear dependence of 
the Drell-Yan production are in the small qx region [§-^|, it is very important to derive a 
quantitative description of the nuclear dependence coefficient a (A, qx) or the ratio R(A, qx) 
for the small qx region. 

Furthermore, J/ip data in hadron- nucleus collisions show that the a (A, qx) (or R(A, qx)) 
for J/ijj production as a function of qx has similar features as that of the a(A, qx) (or 
R(A, qx)) for the Drell-Yan production J7|Q. Recent data from Fermilab experiment E866 
show that a(A,q T ) for J/ip production has an universal shape, but its magnitude depends 
on the range of xf H- J/V suppression in relativistic heavy ion collisions was predicted to 
signal the color deconfinement |I7j . Recently, significant J/ip suppression has been observed 



in existing fixed target experiments []T8| , pL9|| . The data have generated a lot of theoretical 
discussions in searching for the mechanism of the observed J/ip suppression. Understanding 
the features in a(A,q T )j/^p is very valuable for such investigations. 

In this paper, we argue that in the small qx region, a (A, qx) for the Drell-Yan production 
is given by an universal function: a\+hi q T , where a>\ and h\ are completely determined by the 
measured ratio of total cross sections, ^r/^r> an d the averaged transverse momentum 
square (qx)r>Y i n hadron- nucleon collisions, plus perturbatively calculable quantities, such 
as the transverse momentum broadening A{q T )x>Y in hadron-nucleus collisions. The {qx)r>Y 
and A(q T ) DY will be defined in Eqs. (0) and ([12]). 

According to the generalized QCD factorization theorem [|TjJ, like all perturbatively 
calculable hadronic quantities, the transverse momentum broadening for the Drell-Yan pro- 
duction, A(q T )x)Y, can be factorized into a convolution of an infrared safe hard part and 
corresponding universal quark-gluon correlation functions. The infrared safe hard part is cal- 



culable in QCD perturbation theory and was derived in Ref. [20|. Although the quark-gluon 



correlation functions are nonperturbative in nature and unknown, just like the well-known 
parton distributions, these functions are universal and they appear in the factorized expres- 
sions of other physical observables. For example, the same quark-gluon correlation functions 
in the expression for A(q T ) dy appear in the factorized formulas for the a(A, qx) in the large 
qx region [O. Because of the universality of these quark-gluon correlation functions, the 
data on A(q T )x)Y can be used to extract the quark-gluon correlation functions, which can 
then be used to predict a (A, qx) in the large qx region. In this paper, we show that the Fer- 
milab E772 data on A{q T )oY can be used to extract the size of the quark-gluon correlation 
functions, which can then be used to predict the a (A, qx) in the small qx region. 

We further demonstrate that our a(A, qx) for the Drell-Yan production in small qx region 
is naturally connected to the perturbatively calculated a(A, q T ) in the large qx region, and 
we also show that the predicted a(A,q T ) is consistent with the E772 data in both the 
small and large qx regions. Furthermore, we show that although the ratio of the transverse 
momentum spectrum R(A, qx) can have a non-trivial dependence on the atomic weight A, 
the nuclear dependence coefficient a (A, qx) is insensitive to the atomic weight A for normal 
nuclear targets. On the other hand, we argue that in a color deconfined nuclear medium, 
the a(A, qx) becomes very sensitive to the A due to a long range color correlation. 

In addition, we show that a(A,qx) for J/ip production in small qx region is given by 
the same functional form: 02 + &2 q\, and the new parameters 02 and 62 are also completely 
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determined by the the ratio of total cross sections ffm/ffju and the averaged transverse 
momentum square (<?t)j/v, m hadron-nucleon collisions, plus perturbatively calculable quan- 
tities, such as the transverse momentum broadening A(gf,)j/^, in hadron-nucleus collisions. 
Our predictions for OLz/^{A,qr) are consistent with recent data from Fermilab E866. 

The predicting power of our nuclear dependence coefficient function a{A,qr) for the 
Drell-Yan and J/ip production is its universal quadratic dependence on q?, and the fact 
that all parameters are completely fixed by either calculable or independently measurable 
quantities. In addition, we predict that the ct(A, gr) is extremely insensitive to the atomic 
number A for normal nuclear targets while it can be very sensitive to A (or the medium 
size) for a color deconfined medium. 

The rest of this paper is organized as follows. In the next section, we argue that the Drell- 
Yan transverse momentum spectrum in hadron-nucleon and hadron-nucleus collisions should 
be well represented by a Gaussian-like distribution in the small gy region at fixed target 
energies. From such a transverse momentum spectrum, we derive a universal expression 
of a(A, qx) for the Drell-Yan production in the small gr region in Sec. We explicitly 
demonstrate that the ct(A, gr) is insensitive to the atomic weight A, and the a(A, q?) in 
the small qr region is naturally connected to the perturbatively calculated a(A, qx) in the 
large g^ region. We then compare our universal function a(A, q?) with E772 data at all 
qx- In Sec. [TV], we provide a(A, q T ) for J/ip production, and discuss the similarity and 
difference between the a(A, q?) for the Drell-Yan and J/ip production, we show that in a 
color deconfined nuclear medium, the at(A, qx) for J/ip production becomes very sensitive 
to the atomic number A (or the medium size), while it is extremely insensitive to the A for 
normal color confined nuclear targets. Finally, in Sec. [VJ, we summarize our main conclusions, 
and discuss the predicting power of our nuclear dependence coefficient a(A,q T ). 



II. THE DRELL-YAN TRANSVERSE MOMENTUM SPECTRUM 

Depending on physics origins, the Drell-Yan transverse momentum spectrum can be 
divided into three regions, as shown in Fig. |I[ The small qx region corresponds to the region 
where g^ < g^ ~ 1 GeV. The spectrum in this region, labeled by I in Fig. [TJ, is dominated 
by the intrinsic transverse momenta of colliding partons. At leading order in perturbation 
theory, the Drell-Yan transverse momentum distribution is given by 



da hN da hN 



dQ 2 d?q T dQ 2 



8\q T ) . (2) 



The effect of the parton intrinsic transverse momentum can be included by replacing the 
5-function in Eq. (|2]) by its Gaussian-representation, S(q x ) = lim r ^ 75=7 exp[— q 2 j2r 2 \ fl2~If , 
and we obtain 

da(I) =N DY —e-^ 2 (3) 
dQ 2 dq 2 iVDy 2r 2 ' [) 

where the superscript "(I)" indicates the region "I", N^y is a dimensional normalization, 
and r is the width of the Gaussian-like distribution. The physical meaning of the N^y and 
r will be determined later. In principle, A^y and r are non-perturbative quantities. 
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The intermediate qx region, labeled by II in Fig. [l|, corresponds to the region where 
Qx — It < It with = kQ and n ~ 1/3 — 1/2. In this region, the physically observed 
scales: Q and qT are both large enough, and in principle, QCD perturbation theory can be 
used to calculate the short- distance partonic scattering and the Drell-Yan cross section can 



be factorized into the following form [14 



where J2 a ,b sum over & U parton flavors a and b, the \i represents both the renormalization and 
factorization scales, the 0's and ®'s represent the parton distributions and the convolution 
over incoming partons' momentum fractions, respectively. In Eq. (f|), the da ab j dQ 2 dq T is a 
short- distance contribution for partons a and h to produce a Drell-Yan lepton pair, and it 
can be calculated perturbatively in terms of a power series of the strong coupling constant 
a s (n 2 ) with the coefficients that depend logarithmically on the ratio of scales: Q 2 /fi 2 and 
q T /fi 2 . The scale \i can be chosen to be equal to Q or q?, or somewhere between, and the 
factorized Drell-Yan cross section in Eq. (^j) should not be sensitive to a specific choice of 
the scale if the factorized expression is reliable. 

If the Q is much larger than the qx, the perturbatively calculated da ab /dQ 2 dq\ will have 
large logarithms, such as \n(Q 2 /q T ), for any choice of the scale \x. These large logarithms are 
due to the gluon radiation off the incoming partons before the hard collision takes place to 



produce a Drell-Yan lepton pair p2| , |23| . Because of the infrared and collinear singularities 
associated with the massless gluons, we can have two powers of the large logarithms for every 
power of the strong coupling constant a s in the perturbatively calculated da ab '/ 'dQ 2 dq T . 
Clearly, the logarithms, (a s In 2 (Q 2 / q T )) n , can be larger than the unity if the Q is much 
larger than the qx- For example, for Z° production at the Tevatron energies, the ratio 
Q 2 /q T can be as large as 8 x 10 3 for qx ~ 1 GeV and Q ~ 91 GeV. Consequently, the high 
order corrections in a s can be as important as the lower order terms, and therefore, the 
perturbatively calculated da ab /dQ 2 dq T at a fixed order of a" are no longer reliable. 

Using the renormalization group method, the large logarithms, (a s In 2 (Q 2 / q T )) n , in the 
perturbatively calculated da ab / dQ 2 dq\ can be resummed to all orders in a s |2^J2^] . As 
pointed out in Ref . |23| , the resummation of the large logarithms can be systematically car- 



ried out by solving a corresponding renormalization group equation in the 6-space, which is 
the Fourier transform of the transverse momentum g^-space. The kernel of the renormal- 
ization group equation can be calculated order by order in a s in QCD perturbation theory. 
In order to obtain the Drell-Yan transverse momentum distribution, one has to perform 
the Fourier transform of the 6-space solution of the renormalization group equation back to 
the qx space; and it is necessary to have the 6-space solution for all values of b. However, 
since the kernel was calculated in perturbative QCD, the 6-space solution is only reliable 
for small values of b (i.e., b 1/K.qcd)- Therefore, a nonperturbative function F NP (b) 
had to be introduced to cover the large b region | 2~3| . The parameters in the F NP (b) were 
determined by fitting experimental data on the Drell-Yan qx spectrum [pl, p5|1 . The role of 
this nonperturbative function F NP (b) is similar to that of the input parton distributions 
(f)(x, Qq) when we solve the renormalization group equations (or the DGLAP equations) for 
the parton distributions' Q 2 -dependence. 

As argued in Ref. [p3fl , the nonperturbative F NP (b) should have a Gaussian form in the 
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6-space. Since large b values correspond to a small qr region after the Fourier transform, the 
Gaussian form of F NP (b) is consistent with our q? spectrum in region (I). Furthermore, the 
Drell-Yan q? spectrum in region (II) at fixed target energies is mainly determined by the 
resummed 6-space distribution in large b region, which is completely dominated by the shape 
of F NP (b). Consequently, the resummation performed in the small b region does not have 
any noticeable effect on the low qr spectrum at fixed target energies. Experimentally, all 
existing data for the Drell-Yan continuum between J/-0 and T peak at fixed target energies 
HID can be well represented by the Gaussian-like qx distribution for q T < q^. In Fig. ||, we 
show that the Drell-Yan data at both 800 GeV |2B| and 400 GeV |27J can be well represented 
by a Gaussian-like fit for q T up to 2.5 GeV, and Q = 5.5, 8.5, and 11.5 GeV, respectively. 
The widths of the Gaussian-like fits in Fig. are consistent with the parameters in the 
F NP (b) in the resummed Drell-Yan spectra p4| , p5f . 

Actually, this finding should not be too surprising because for a typical Drell-Yan pair 
of Q ~ 5 GeV and g T ~ 1 GeV at fixed target energies, the ratio Q 2 jq\ is about 25, which 
is much smaller than the 8 x 10 3 for a Z° production at the collider energies. On the other 
hand, when Q 2 is much larger than the typical gf-, the gr spectrum in this region has an 
enhanced peak at a finite value of qr and a slower fall off in comparison with the Gaussian- 
like distribution, as sketched in Fig. [TJ. The size of the enhancement and the size of the region 



depend on the value of Q and the range of qx [25]. Furthermore, because of the slower fall 
of the transverse momentum distribution with a full resummation of the large logarithms 
at low qr, the transition between the resummed distribution and the large qr distribution 
calculated in the fixed order perturbative QCD should be smoother than the unphysical 
kink, as shown in Fig. |2], at the point where the Gaussian-like distributions are matched 
with the large qx perturbative tails. Since we concentrate on Drell-Yan q? spectrum at fixed 
target energies in this paper, we will effectively merge regions (I) and (II) in the rest of our 
discussions. 

The large qr region, labeled by III in Fig. [TJ, corresponds to the region where qx > qj>- 
In this region, the transverse momentum spectrum calculated in conventional fixed order 
perturbative QCD should be reliable [p8 , p9 |, as shown by the solid line in Fig. |l]. This 



perturbative tail of g^-spectrum has a typical power-like behavior. 

In conclusion, the Drell-Yan g^-spectrum at fixed target energies can be represented by 
a Gaussian-like distribution in the small qr region and a perturbatively calculated tail in 
the large qr region, 

- <tr) + ItL l %r - Qt) 



dQ 2 dq 2 dQ 2 dq 2 w ^' dQ 2 dq 2 



dQ 2 dq\ \ dQ 2 dq\ dQ 2 dq, 



2 



9(q T -qL), (5) 



where is the perturbatively calculated gr-spectrum p8| , |29[| , which is valid for region 

(III), and J^dq-z is t ne Gaussian-like distribution defined in Eq. (^|), which is able to fit data 

in regions (I) + (II). Since ^2^1 is calculable, q^ and parameters Ndy and r in Eq. @ are 
the only unknowns for the Drell-Yan g^-spectrum at fixed target energies. 

Taking the moments of the Drell-Yan q^ spectrum in Eq. (|5]), we can express the pa- 
rameters Ndy and r in Eq. @ in terms of physically measurable quantities. For example, 
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by integrating the qx spectrum in Eq. fl5|), we obtain the Drell-Yan total cross section as 
the zeroth moment of the spectrum. Using the data shown in Fig. ||], we found that the 

contribution from the second term in Eq. (|5|) J q L {jt^zj^r — Jq^^ ) ^Qt * s muc h less than 

one percent of the contribution from the first term / f ^§^2 J dq\. Therefore, up to the less 
than one percent uncertainty, we obtain the normalization Noy in Eq. (W) as 



da 



>DY 



dQ 2 

Define the Drell-Yan averaged transverse momentum square as 



(6) 



/ 2 \ _ [ ' ql{d(j/dQ 2 dql)dql 

{(!t) - daJdQ 2 • [V 

Substituting Eqs. (|) and (|6|) into our definition in Eq. ([/]), we obtain 

(q 2 T ) = 2r 2 + T(q^T 2 ) , (8) 



with 



da/dQ 2 J ^ * J \dQ 2 dq 2 , dQ 2 dq. 

In principle, T(q^,r 2 ) is calculable, and its r 2 dependence is from the da^ 1 ' j l dQ 2 dq\. Even 
though (gf,) enhances the contributions from the perturbative tail, we find that the T(q^, r 2 ) 
in Eq. @ is much less than ten percent of the first term for the data in Fig. |2] at fixed target 
energies. Therefore, by iteration, we derive the width of the Gaussian-like distribution 
da^jdQ 2 dq\ in Eq. (g) as 

2r 2 ^(g 2 )-r(g|,r 2 = (g 2 )/2), (10) 

where the precise value of q^ will be determined later. 

Substituting Eqs. (|5]) and (^) into Eq. (|3|), we obtain the Drell-Yan g T -spectrum at small 
qx in hadron-nucleon collisions as 



dQ 2 dq 2 dQ 2 (q 2 ) h D N Y - T(q^ ' 1 ) 



where the superscript hN represents the hadron-nucleon collisions, (qfyoY * s the averaged 
transverse momentum square for the Drell-Yan production in hadron-nucleon collisions, 
and the short-handed notation T(q^) h D T Y = T(q^r 2 = {q 2 ) h D T Y /2) h D ^ with T(q^T 2 ) defined 
in Eq. (|9|). We emphasize that all quantities, da hN /dQ 2 , (q^DY, and r(g^)^-, which 
completely define the Drell-Yan g^-spectrum in the small qr region, are either independently 
measurable or perturbatively calculable. For the Drell-Yan data in Fig. ^ at fixed target 
energies, we found that (q^DY ranges from 1.5 GeV 2 to 1.8 GeV 2 , and r(g^)^ is less or 
about ten percent depending on the value of q^. 

In order to derive an expression for the nuclear dependence coefficient a (A, q?), we need 
to know the Drell-Yan transverse momentum spectrum in hadron-nucleus collisions. In 
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hadron-nucleus collision, initial-state partons have much larger probability to have multiple 
scattering before the hard collision to produce the Drell-Yan pair. Therefore, in order to 
obtain the transverse momentum spectrum in hadron-nucleus collisions, we need to consider 
contributions from both the single scattering and the multiple scattering. For the single 
scattering contribution to the transverse momentum spectrum in hadron-nucleus collisions, 
we can use the same arguments given above for the spectrum in hadron-nucleon collisions. 
Therefore, at fixed target energies, the single scattering contribution to do hA jdQ 2 dq\ in the 
small qx region can be well represented by a Gaussian-like distribution. 

For the Drell-Yan production, parton initial-state interactions dominate the multiple 
scattering because the virtual photon does not interact strongly. The initial-state multiple 
scattering can provide extra transverse momentum to a parton participating the hard col- 
lision, and consequently change the transverse momentum spectrum of the produced Drell- 
Yan pair. Therefore, multiple scattering leads to a broadening in the averaged transverse 
momentum square (qj) hA , which is defined in Eq. (0) with the superscript hA representing 
the hadron-nucleus collisions. Define 

A<?£> = (q 2 T ) hA ~ (q 2 T ) hN , (12) 

and A(g|>) is often called the nuclear broadening of the averaged transverse momentum 
square. It was demonstrated in Refs. p0| , |30| that the nuclear broadening A(g|>) is calculable 



within QCD perturbation theory. In terms of the double scattering and the generalized 
factorization theorem []15| |, the transverse momentum broadening for the Drell-Yan pairs, 



^■(Qt)dy, was calculated at the leading order in a s , and was expressed in terms of the 
quark-gluon correlation functions inside a nucleus IpCf] : 

2 = (4nW\ E q e 2 q Jdx'<P, /h (x')T q f(Q 2 /sx')/x' 
[(lT)DY { 3 )' E.e^Jdx'^x^^/AiQysx^/x" [ } 

where J2q runs over all quark and antiquark flavors, e q is the fractional charge of quark and 
antiquark of flavor q, and a s is the strong coupling constant. In Eq. (HJ), the s = (p + p') 2 



with p' the four momentum of the incident hadron and p the averaged momentum per 
nucleon for the nucleus, and 4>q/h and (p q /A are parton distributions of the hadron and the 
nucleus, respectively. The T^ g H (x) of flavor q in Eq. (|13|) is the "soft-hard" quark-gluon 
correlation function in a nucleus and is defined by [ |10| , [20|j 

x 1 (P A \F a + (y^ q (0)l + My-)F +a (yl)\PA) , (14) 

where we assumed that the nucleus is moving along the +z direction, and the superscripts 
"+" and "— " represent the normal + and — components in a light-cone coordinate, respec- 
tively. In Eq. (|T4j), the y, y\ and y2 represent the position coordinates of the quark and 
gluon field operators. Although these quark-gluon correlation functions are nonperturbative 
and unknown, they are universal and as fundamental as the well-known parton distributions 
||16|| . They appear not only in Eq. (|1^) for the nuclear broadening A(g|,)oy, but also in the 
nuclear dependence of the Drell-Yan g^-spectrum in hadron-nucleus collisions in the large qx 



S 



region |Tl|,|T^]. In principle, by using the nuclear dependence of the Drell-Yan g^-spectrum 
in the large qr region, we can extract these quark-gluon correlation functions and predict 
the nuclear broadening A^^dy, or vise versa. 

On the other hand, by comparing the operator definitions of the quark-gluon correlation 
functions T^ g (x), defined in Eq. fll4]), with the definitions of normalized quark distributions, 



1 A (x) = -(j) q/A {x) 



1 

A 



dy 

2tt 



(15) 



one finds that the only difference between the T^ H (x) and q (x) is an operator factor 



dyi dy 2 
2tt 



0{y--vI)0{-ih)F + a{ih)F +a {vI 



If these two gluon field strength F's are close together in a color singlet nucleon in a normal 
color confined nuclear target, which limits the / dy± (or / dy 2 ) to a nucleon size, the extra 
integration / dy 2 (or / dy±) can be extended to the nuclear size to give an extra A 1 / 3 factor 
to the Tq g H (x) over q A (x). Therefore, one can estimate the size of the quark-gluon correlation 
functions by modeling them as [H^fTo] : 



T. 



SH, 
19 ' 



(16) 



where A is a universal parameter proportional to the averaged color field strength square 
in position space in a normal nuclear target, and the q (x) are the effective nuclear quark 
distributions normalized by the atomic weig ht A. With this model for T^ H , the Drell-Yan 
transverse momentum broadening A^^dy can be expressed as [20 



A(ql) DY 



A 2 A 1/s = b DY A l/3 



(17) 



The factor of A 1 / 3 in Eq. 
nuclear matter 
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[17|) is an unique feature of parton double scattering in a normal 
Data for the transverse momentum broadening is consistent with the 



features of double scattering @J20 



Same as the single scattering case, double scattering contributions to the Drell-Yan qx- 
spectrum in the small qx region is also dominated by the intrinsic transverse momenta of 
scattering partons. Similar to Eq. (Q), at the leading order in perturbation theory, the double 
scattering contributions to the Drell-Yan da hA jdQ 2 dq\ can be described by fl20| 



da h D A 
dQ 2 d 2 qr 



oc T qg (x, X!,x 2 , k T ) S 2 (q T - k T ) 



(18) 



where the subscript D indicates the double scattering, and the proportionality coefficient 
depends on the structure of partonic double scattering. In Eq. (0), qr is the transverse 
momentum vector of the observed Drell-Yan pair, represents the intrinsic transverse 
momentum vector of the gluon which gives additional scattering, and T qg (x, x±, x 2 , hr) rep- 
resents a general quark-gluon correlation function with x, x±, and x 2 being the light-cone 
momentum fractions carried by the quark and gluon fields, and fcy is the magnitude of kx 
16| , |20[| . Integrations over parton momenta dxdx\ dx 2 reduce the general T qg (x, xi,x 2 , kr) to 
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either , defined in Eq. ([14]), or T® H which is defined later in Eq. (p5|) , depending on the 
structure of the partonic part of the double scattering. Similar to Eq. (||), we can broaden 
the 5-function in Eq. (|18|) by its Gaussian-representation to include the effect of intrinsic 
transverse momentum of incoming quark or antiquark. We can also take into account the 
effect of gluon's intrinsic transverse momentum by using a Gaussian-like distribution for 
the fcr dependence of the T qg . After integrating over the intrinsic fey-dependence, we find 
that just like the single scattering case, the double scattering contributions to the Drell- 
Yan g^-spectrum in the small region can also be expressed in terms of a Gaussian-like 
distribution defined in Eq. with a broadened width. 

In order to understand the Drell-Yan spectrum in the region (II) in hadron-nucleus 
collisions, in principle, we have to carry out the resummation of the large logarithms at 
the presence of multiple scattering. Following the generalized factorization theorem ||15|| , 
the resummation of the large logarithms at the presence of factorizable partonic multiple 
scattering can be carried out by using renormalization group method. In order to resum the 
large logarithms due to gluon radiations, it is necessary to calculate anomalous dimensions 
of the multiparton correlation functions for deriving the kernel of the renormalization group 
equation. However, following the same arguments given above for the single scattering 
case, the resummed Drell-Yan g^-spectrum in region (II) at the fixed target energies is 
mainly determined by the nonperturbative input distribution F NP (b), which generalizes 
the perturbatively resummed Drell-Yan 6-space distribution to the nonperturbative large b 
regime in order to perform the Fourier transform back to the gr-space. As argued above, the 
Drell-Yan qx spectrum at the presence of multiple scattering should still have a Gaussian-like 
distribution in the small qx region, and we expect that the nonperturbative F NP (b) at the 
presence of multiple scattering should also have a Gaussian-like distribution. Therefore, at 
the fixed target energies, we can merge the region (I) and (II) of the Drell-Yan g^-spectrum 
in hadron-nucleus collisions. Since we limit this paper to test the consistency of QCD 
treatment of multiple scattering at fixed target energies, which has not been done, we will 
defer a complete analysis of QCD resummation at the presence of multiple scattering to a 
future publication, which will be very important for the nuclear dependence of Drell-Yan q? 
spectrum at the collider energies. 

Similar to the case in hadron-nucleon scatterings, when q^ > g^, the Drell-Yan q^- 
spectrum should have a perturbatively calculable tail in region (III), and its nuclear de- 
pendences are perturbatively calculable and were calculated in Ref. [fyj. Therefore, the 
Drell-Yan g^-spectrum in hadron-nucleus collisions should have the same form as that in 
Eq. (|5|). Using the calculated perturbative tail for the Drell-Yan spectrum in hadron-nucleus 
collisions, we find that similar to the hadron-nucleon case, the second term in Eq. (|5|) con- 
tributes only about one percent of the total cross section da hA /dQ 2 . Therefore, similar to 
Eq. (H), we can identify the normalization factor of the Gaussian-like distribution to be the 
total cross section da hA /dQ 2 . 

Similar to Eq. flTTD, we derive the Drell-Yan transverse momentum spectrum at small qx 
in hadron-nucleus collisions as 

da hA da hA 1 -&/««3.>&-r(g£)&) (19) 



dQ 2 dq\ dQ 2 {q¥) h A-n^) h DY 

where the superscript hA represents the hadron-nucleus collisions, (q^DY = (Qt)dy + 
A(q^) DY , and the nuclear broadening A(q^) DY is given in Eq. (pf) or Eq. (|17|). In Eq. (|19|), 



10 



r(gf:)^ = r(g£,r 2 = (gf)M7 2 )i)y with v ^ q L^hA defined in Eq . (|) with the ca l cu l a _ 
ble hadron-nucleon g^-spectrum replaced by the corresponding hadron-nucleus g^-spectrum. 
Similar to the hadron-nucleon case, we find that the r(g^)^ i n Eq. ( |19D is much smaller 
than corresponding (q T )rjY- Since the Drell-Yan total cross section da hA /dQ 2 has almost no 
nuclear dependence, Eq. ( |i9|) shows that the nuclear broadening A{q T )c Y determines the 
Drell-Yan g^-spectrum at small qx in hadron-nucleus collisions. 

From Eqs. (|TT]) and (|T^), we can derive the nuclear dependence coefficient a(A,qT) for 
the Drell-Yan process in small qx region. 



III. THE NUCLEAR DEPENDENCE COEFFICIENT FOR THE DRELL-YAN 

PRODUCTION 

In this section, we derive the nuclear dependence coefficient a(A, qx) for the Drell-Yan 
production in the small qx region. We also show that the derived a(A, qx) in the small qx 
region is naturally connected to the perturbatively calculated a (A, qx) in the large qx region. 

Substituting Eqs. (|ll]) and (^) into Eq. ([!]), we derive the aDY(A,qx) for the Drell-Yan 
production in the small qx region: 



a DY (A,q T ) = 1 + 



1 



ln(A) 



In (R DY (Q 2 )) + In 



1 



1 + Xdy, 



Xdy 



q\ 



1 + Xdy 



J2\hN r(„L\hN 

It) dy 1 yqrJDY 



where R% Y (Q 2 ) = (l/A)(da hA /dQ 2 )/(da hN /dQ 2 ). The X dy in Eq. is defined by 



Xdy = 



A(q 2 ) 



DY 



AT{q^) DY 



A(q 2 } 



DY 



(4)^ - m) D N Y 



TIDY 



(20) 



(21) 



with A{q T ) DY given in Eq. flT7|), and AT(q^) DY = ^(q^DY ~ T(Qt)dy- Because T(q^) DY 
contributes less than ten percent to {qx) dy in both hadron-nucleon and hadron-nucleus 
collisions, AT faff) dy is also much smaller than A(q T )x>Y- Since the xdy and (gf-)zjy have 

01) for the 



no dependence on qx, the nuclear dependence coefficient (Xdy(A, qx) in Eq. 
Drell-Yan production in the small qx region shows an universal quadratic dependence on 
qx- axjY(A,q T ) = a\ + biq T . The parameters ai and b\ are completely determined by 
either perturbatively calculable quantities (such as A(q T ) and the T(g|;)'s) or independently 
measurable quantities (such as R^ Y (Q 2 ) an d (qx) hN )- 

Using E772 [g and NA10 @ data on A(q T ) DY and Eq. (0), we obtain that A(q T ) DY ~ 
0.022V4 1 / 3 GeV 2 , which is between 0.05 and 0.13 GeV 2 for all practical nuclear targets. On 
the other hand, the value of (q T ) hN for the data in Fig. |2] is ranging from 1.5 to 1.8 GeV 2 . 
Therefore, the xdy defined in Eq. ( |2"TD should be a small number. 

~~ ]), and taking the small xdy limit, we obtain 



Substituting Eqs. (|T7|) and ( pT|) into Eq. 



a DY (A,q T ) wl + 



1 + 



Xdy 
ln(A) 

bpY 

( n 2\hN 
yIti DY 



-1 + 



q 2 r 



( n 2\hN _ 
WT7 DY 

q\ 

" ln 2 \ hN 
yIT! DY 



T(q, 



TIDY 



+ o( 
KM*)) 



(22) 
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where buy = (47r 2 a s /3)A 2 is defined in Eq. (|T7D. Here we used the experimental fact that 
Rdy(Q 2 ) ~ 1- in deriving the second line in Eq. (p2j) , we used A 1 / 3 ~ ln(A), which is a 
good approximation for most relevant targets, and we neglected the small ^(<1t)dy- F° r a 
given beam energy, (gf^y f° r the Drell-Yan pairs has a very week dependence on Q and no 



dependence on A. Therefore, Eq. (p2[) shows that the leading contribution to the any (A, qr) 
for the Drell-Yan production is a function of qr and independent of A in the small qr region. 

From the size of the correction term in Eq. ([2"2"D, it is clear that the acy(i,gr) for the 
Drell-Yan production is insensitive to the atomic weight A of targets. Furthermore, from 
Eq. (JZUp, we obtain 



aDY{qr = 0) « 1 



ln(l + xpy) 
In(A) 



ln(A) 



1 - 



(4) h D N Y 



(23) 



From Eqs. ( gl]) and fl23|) , it is clear that C(dy((1t = 0) is determined by the nuclear broadening 
A(<7t)dk- If there is nuclear broadening in averaged transverse momentum square (i.e., 
Xdy > 0), the transverse momentum spectrum is suppressed at qx ~ 0. For a normal 
nuclear target, the amount of the suppression at qx ~ is insensitive to the atomic number 
A of targets. 

Eq. fl2"0|) shows that as qx increases, ooy(A qx) increases as well, and it changes from 
less than one to a value larger than one. Let q T correspond to the critical value of gr at 
which a£>y(A Qt) = 1' anc ^ we h &ve 



{q T f « (g 



TIDY 



1 + Xpy 



ln(l + XDy ) « (g 2 )^ [1+OM] 



(24) 



Like «£)y(A, gr = 0), the q T in Eq. (^) has a very weak dependence on the atomic number 
A for a normal nuclear target. 

In the rest of this section, we show that the a>DY{A,qT) in Eq. ( f2~0D (or approximately 
in Eq. (E^D) for the small gr region is naturally connected to the perturbatively calculated 



ar,Y{A^q T ) in the large qx region JTT 



The nuclear dependence of the Drell-Yan transverse momentum spectrum in the large 
qx region can be calculated in QCD perturbation theory |11|,12|. The ratio of the Drell-Yan 

A a(A, qT )-i defined in Eq. (0) was derived in Ref. for the large 
that at the leading order in a s , the R(A, q T ) depends on two 
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g^-spectrum R(A, qx) = 
qr region. It was found 
types of parton correlation functions due to different structure of partonic double scattering. 
For example, consider a double scattering Feynman diagram in Fig. |3], which contributes to 
the nuclear dependence of the Drell-Yan g^-spectrum in the large gr region. In Fig. ^, x, 
Xi, and X2 are independent parton momentum fractions for the double scattering. For the 
contour integration of dx^s, taking the residues of poles of the partonic double scattering 
diagram can eliminate the oscillation factor e~ XiP+y i for corresponding integrations dy~^s in 
position space [|l(],|ll|]. The lack of the oscillation factor for integration of position variable 
?/~'s leads to an A 1//3 -type (or nuclear size) dependence. The partonic double scattering 
diagram in Fig. |3] has two types of poles: one type that corresponds to x — x\ ^0 and 
x — X2 — ► 0, and the other type that corresponds to X\ = x 2 = x a and x — x\ = x — x 2 = Xb 
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with both x a and x b finite. The first type poles lead to the so-called "soft-hard" quark-gluon 
correlation functions defined in Eq. fll4j) , which have both dyi and dy 2 integrations free of 
the oscillation factor. The second type poles lead to the so-called "double-hard" quark-gluon 



correlation functions 10,11 



T g D g H (x a ,x b ) = J-f^Lf^ir p+dy 2 e(y--y;)9(-y 2 )e^ + y-e^ + ^-y^ 



qg v XbP + J 2n J Air « 

x (P A \F a + ( y2 -)$ q (0h + i; q (y-)F + "(y-)\P A ) , (25) 

where all the variables and the superscripts carry the same meaning as those in Eq. ([Tip. 
Unlike the soft-hard correlation functions, only one dy 2 (or cfa/f) integration in the double- 
hard correlation functions defined in Eq. (|25| ) is free of the oscillation factor. But, this dy 2 
(or dy^) integration is just enough to produce an A 1//3 -type enhancement to the Drell-Yan 
gy-spectrum pH|JTT[ ]. In addition to the double-hard quark-gluon correlation functions, the 
nuclear dependence at large qx also receives contributions from quark-quark and gluon-gluon 



correlation functions ||11|| . By comparing the operator definitions of the normal quark and 
gluon distributions with the operator definitions for the double-hard correlation functions, 
similar to Eq. (|16|) , the double-hard quark-gluon correlation functions can be parameterized 



as 31 



T q D g H (x a , x b ) = (2nC) v4 4/3 q A {x a ) G A (x b ) , (26) 

with q (x a ) and G A (x b ) the effective nuclear parton distributions normalized by the atomic 
number A. The constant parameter C represents the size of quantum interference between 
different nucleon states and possible overlap of two nucleon parton distributions in a nu- 
cleus. Currently, there is no direct observable yet to extract information on T DH [12] or to 
determine the parameter C in Eq. (|26|) . 



When x\ (or x-i) goes to zero, the corresponding parton fields reach the saturation region, 
and the T DH is reduced to T SH . Therefore, if we assume the full quantum interference (in 
saturation region), the proportionality constant C in Eq. ( ]2T)| ) can be related to the A 2 in 



Eq. ( fig) by using the following identity: 



lim x b T q D g H (x a , x b ) = T q f(x a ) , (27) 



x b — >Q 

From Eq. (p7|) and the definitions of T® H and T q H , we obtain 



A 2 

C ~2nxG(x)U > (28) 
where xG(x)\ x ~o is the gluon momentum density at x — > limit, which is believed to be 



order of one at low energy p2| , p3| . The precise value of xG(x)\ x &o depends on the physics of 
parton saturation [Q, and it has not been well measured experimentally. In this paper, we 
choose xG(x)\ x ~o ~ 3 based on the simple form G(x) ~ 3x _1 (1 — x) 5 , which is a result of 
a simple power counting for gluons at large x, Regge behavior for gluons at small x, and an 
assumption that gluons carry 50 percent momentum of the beam hadron. Different choice 
for xG(x)\ x ~o results into a different overall normalization factor C. The A 2 in Eq. (p8[) is 



the same as the one in Eq. fll6|) for the soft-hard quark-gluon correlation functions. Other 
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correlation functions can also be determined in a similar way. Therefore, the A 2 becomes 
the only parameter to fix the normalization of both the soft-hard T SH and the double-hard 
T DH correlation functions. If we extract the value of A 2 from the data on the transverse 
momentum broadening A(gf,) DY by using Eq. fllTf), we can fix our predictions of the nuclear 
dependence for the Drell-Yan g^-spectrum in the large qr region. 

On the other hand, by assuming that the quantum interference between parton field 
operators from different nucleon states are strongly suppressed, one derives 

C = 0.35/(87rr 2 ) GeV 2 (29) 

with ro ~ 1.1 — 1.25 fm, which is just a geometric factor for finding two nucleons at the 



same impact parameter [11,31 



Because of a combination of a small value of the measured A 2 from Drell-Yan data |||4j 
and a choice of [xG(x)] x ~q ~ 3, the numerical values of C determined by Eqs. ( |2"9"D and ( |2"8"| ) 
approximately differ by a factor of 20. We believe that the numerical value for the C given in 
Eq. fl2"9"|) without any quantum interference between parton field operators should represent a 
possible maximum for the C, while the value given in Eq. (|28"D with full quantum interference 
(in saturation region) represents a possible minimum for the C. The parameter C reflects 
the size of quantum interference involved in the multiparton correlation functions, and its 
precise value should be between those obtained by Eqs. (|29"D and (|^). More discussions on 
the consequences of a different value of the C are given in Sec. 0. 

In Fig. f|, we plot the arjY(A,q T ) as a function of qT for various nuclear targets. For 
the small qx region, we used our ajjyiA, q T ) in Eq. (p0|). For the large qx region, we used 
the perturbative calculations of R(A,qx) = A a ^ A ' qT ^~ 1 in Ref. |TI]. In plotting Fig. |]a, we 



used {q T ) h ^ Y = 1.8 GeV 2 , and A(q T ) DY = 0.022A 1 / 3 GeV 2 . We choose the value of A 2 to be 
consistent with boy = 0.022 GeV 2 . The solid and dashed lines correspond to the nuclear 
target C (A=12) and W (A=184), respectively. For the large qr region, we plotted the 
perturbatively calculated anyiqr) [0 with the maximum and minimum values of the C 



discussed above. It is clear from Fig. (|a that our nuclear dependence coefficient for the 
Drell-Yan g^-spectrum in the small qx region is naturally connected to the perturbatively 
calculated aoY(A,qT) in the large qx region. Without any adjustment of parameters, we 
choose the q^ to be the transverse momentum at the cross point between the low qT and 
high qT spectrum. Without any extra free fitting parameters, our predictions shown in 
Fig. ^ are consistent with data in small q? region, and due to the large error in data at 
high qT, current Drell-Yan data at the fixed target energies are consistent with almost any 
value of the C between the maximum and the minimum discussed above. Different value of 
C resulted in different value of q^ at the crossing point. From Fig. |]a, the value of q? is 
between 2.3 — 4.0 GeV, which is very reasonable. 

As expected from Eq. (p2]), the a D y(A,g T ) in Fig. [|a shows very little dependence on 
the atomic weight A. Therefore, the ocdy{A,Qt) for the Drell-Yan production in the small 
qT region has scaling on A and shows an universal quadratic dependence on q?. The scaling 
violation due to the potential A-dependence is extremely small for a normal color confined 
nuclear target. Since (q^DY can depend on the beam energy and Q as well as the scaled 
longitudinal momentum (xp), our universal function aTjy(A, qx) in Eq. (|20"D can depend on 
these observables accordingly P0|. 



In Fig. f|b, we compare our otDY{A,qT) shown in Fig. [|a with data from Fermilab ex- 
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periment E772 PJ3^]. Using the aDY(A,qx) in Fig. f|a, we plot the ratio of cross sections 
R(A, qx) = A a ( A ' 9T )~ 1 as a function of qx- All parameters are the same as those used in 
plotting Fig. [|a. Also plotted in Fig. |]b are data from E772 on four different nuclear tar- 
gets: C, Ca, Fe, and W [|,@]- It is very clear that without any extra fitting parameter, our 
predictions for the nuclear dependence coefficient C(dy(A, qx) provide a good description of 
the data on R(A, qx) for all values of qx- The A-dependence in R(A, qx) in Fig. |b is clearly 
consistent with the ^-independence of (Xdy{A, qx) shown in Fig. |]a. Better data can provide 
even more critical tests of our predictions for ai DY (A, qx)- 

Fig. ^ demonstrates a clear consistency of the generalized QCD factorization theorem |T5| . 
which was used to calculate the nuclear broadening A(q T ) DY and the nuclear dependence 
in the large qx region. Even though the nuclear broadening A(q T )x)Y is dominated by 
information at low qx while the perturbative tail (or the aoY^A, qx) in the large qx region) is 
controlled by dynamics at large qx, our derived axjY^A, q T ) for the small qx region is naturally 
linked to that for the large qx region without free fitting parameters. More discussions on 
the uncertainties, such as the resummation and the choice of the parameter C, will be given 
in Sec. |V|. 



IV. THE NUCLEAR DEPENDENCE COEFFICIENT FOR 3 /if) PRODUCTION 

A major difference between the Drell-Yan production and J/ip production is the final- 
state interaction. Because of the virtual photon, the Drell-Yan production has only initial- 
state interactions while J/ip production has both initial-state and final-state interactions. 
Therefore, one may conclude that the nuclear dependence of the transverse momentum 
spectrum for J/-0 and the Drell-Yan production could be different due to the strong final- 
state interactions. However, we argue in this section that in the small qx region at fixed target 
energies, J/ip transverse momentum spectrum can also be well represented by a Gaussian- 
like distribution, and the nuclear dependence coefficient a(A, qx) for J/0 production is very 
similar to that for the Drell-Yan production. 

Kinematically, hadronic J/ip production is very much like the Drell-Yan production with 
Q ~ Mjm. Therefore, the J/ip transverse momentum spectrum also has three similar 
characteristic regions shown in Fig. |]. However, because J/0 mass Mj/^, is smaller than any 
typical Q measured from the Drell-Yan continuum, the logarithm (a s \n 2 (Mj/^/q T )) n for 
J/i/j production should be less important than any Drell-Yan production. We then expect 
that at fixed target energies, a Gaussian-like distribution can fit J/-0 transverse momentum 
spectrum at small qx even better. In Fig. [|, we plot data on J/t|i transverse momentum 
spectrum from Fermilab experiment E771 [35|. The solid line in Fig. [5] is a fit using a 
Gaussian-like distribution 

^|U w _Le-^, (30, 

where superscript "(G)" stands for a Gaussian-like distribution, Nj/y and r are the normal- 
ization and the width of the Gaussian-like fit, respectively. As expected, the x 2 fo r fitting 
the J/-0 spectrum in Fig. |5] is about 0.9 per degree of freedom, which is much smaller than 
the x 2 f° r fitting Drell-Yan data in Fig. [2|, where the x 2 P er degree of freedom is ranging 
from 1.0 to 2.0. 
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Similar to Eqs. (0) and (0), at the lowest order in perturbation theory, both the initial- 
state and final-state double scattering contribution to the J/ip transverse momentum spec- 
trum in hadron-nucleus collisions is proportional a ^-function: 5 2 (qr — k?), where is the 
total transverse momentum of any produced pre-J/V> partonic state (e.g., a state of charm 
and anticharm quark pair), and kr represents the intrinsic momentum of the gluon which 
gives additional scattering. Following the same arguments that follows Eq. fllq) in Sec. |1[ we 
conclude that in the small qx region, J /if) transverse momentum spectrum in hadron-nucleus 
collisions can be well represented by a Gaussian-like distribution shown in Eq. (0), and the 
normalization iVjM and the width r are given by 

2r 2 «(^-r(^, (31) 

where is the J/0 total cross section, and ((Zt)ja/, is the averaged transverse momentum 
square, defined in the same way as that for the Drell-Yan (q^DY in Eq. (^). In Eq. (j3ll) 

\hN 



F{Qt)j%> is perturbatively calculable and defined as 



,1AM _ 1 /" 2 / dC7 J/^ ^J/V* | .2 



where da^/dq T is a perturbatively calculable g^-spectrum at large q?, and dafj^/dq T is 
defined in Eq. (^) with 2r 2 = (q^jf^- The corresponding transverse momentum spectrum 
in the small region is given by 



da hA 



dq\ aj/ U4)^-m)if/ ■ { } 

The r(g^)j^y, is the small contribution to {q T )jf^ due to the difference between the pertur- 
bative tail and an extended Gaussian-like distribution, and is defined in the same way as 
that in the Drell-Yan case in Eq. (0). 

Introduce the nuclear broadening of J/0 transverse momentum spectrum A(q T )j/^, 

where (qfyjM is the averaged transverse momentum square of J/0 production in hadron- 
nucleon collisions. Similar to Eq. (0), we derive the nuclear dependence coefficient a(A, qr) 
for J /if) production in the small qx region, 



aj/^(A, q T ) = 1 + 



\n(A) 



In (Rf /ip ) + In 



1 + Xj/V, 



^ XJ/i> Qt 



>j/-<p wr;j/t/i 

MA I ^hN 



(35) 



where Rf /ip = (l/A)a^/a^, and is defined by 
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fixed by either perturbatively calculable quantities (such as A(gf,)j/^ and r(g^)j/^'s) or inde- 
pendently measurable quantities (such as the ratio of the total cross section (1/A)a^f^/ a^ N 



= A(gj,) J/v ,- Ar(gjj)j/^ _ Hips/* 

AJ/-0 — l„2\hN r(„L\hN ~ /„2\WV ' V dD / 

WT/J/V> — 1 WTjj/V WT/J/V> 

Similar to the Drell-Yan case, r(gr)j/^ and Ar(g^)j/^ are perturbatively calculable, and 
much smaller than (g|.}ju and A(gf.)j/^, respectively By comparing Eq. (|20| ) with Eq. (p5|), 
we can see that the oij/^A, gr) for J/?/> production has the same universal quadratic depen- 
dence on qT'. aj/^(A,qT) = a<i + b 2 q^, and the parameters a 2 and b 2 are again completely 

#) Jy ty andr(g£)j/v>' 

"'J/V/ U J/V 

and the averaged transverse momentum square (qfyjhp)- 

For hadron-nucleus collisions at fixed target energies, clear nuclear suppression for J/ip 
total cross sections has been observed [T^|T^|. Various theoretical explanations of the ob- 
served J/ip suppression have been proposed |5fJ. Unlike the virtual photon in the Drell-Yan 
production, the produced cc are likely to interact with the nuclear medium before they exit. 
If one assumes each interaction between the cc pair and the nuclear medium is about the 
same and can be treated independently, one naturally derives the Glauber formula for the 
suppression 

where (3 is a parameter characterizing the size of suppression, and La oc A 1 ^ 3 is the effective 
length of the nuclear medium. In terms of Glauber theory of multiple scattering, the j3 = 
PN^abs with the nuclear density pn and an effective absorption cross section <7 a b s for breaking 
a J/ip meson or for changing a coherent pre- resonance cc pair to a pair of open charms. The 
characteristic feature of the Glauber formula is a straight line on a semi-log plot of (Tj^, vs. 
La- With a large effective <r a b s ~ 6 — 7 mb, the Glauber formula given in Eq. ( p7[) fits almost 
all J/V suppression data in hadron-nucleus and nucleus-nucleus collisions f36j, except the 



strong suppression observed in Pb-Pb collisions ||19|| . 

If one believes that there is no color-deconfinement or quark-gluon plasma was formed at 
present heavy ion collisions at fixed target energies, one has to overcome at least two obvious 
difficulties of the Glauber formula for the suppression: the size of the effective absorption 
cross section cr a b s , and the observed non-linearity of the suppression on the semi- log plot of 
vs. La- It was argued in Ref. ]37| that because of the size of the cc pair produced in 
a hard collision, the J/ip meson has to be formed several fermis later after the production 
of the cc pair. Therefore, the suppression of J/tp in hadronic collisions should be a result 
of multiple scattering between the produced cc pairs and the nuclear medium before the 



pairs exit medium [38,39] . Recall the fact that the Glauber formula in Eq. (|37j) is derived 
for the propagation of a single particle in a medium. On the other hand, the pre-resonance 
J/i/j is a state of at least two particles (e.g., a cc pair), and such a two particle state can 
be changed when it goes through the nuclear medium due to the multiple scattering. Just 
like the random walk for a pair of particles moving through a medium, multiple scattering 
increases the relative momentum between the pairs, and therefore, decreases the effective 
phase space for the pairs to form a bound state J/ip meson. The observed suppression should 
be an immediate consequence of the fact that a cc pair with a larger relative momentum is 
more likely to become a pair of open charm mesons than a bounded J/ip. In this picture of 



17 



3 ft/) suppression, the absorption cross section cr a b s in Eq. (]37[) is an effective cross section for 
breaking a cc pair of relative momentum Ak c5 into a pair of open charm mesons, and clearly, 
such an effective cross section depends on the value of Ak c5 . The larger Ak c5 is, the larger 
the effective cr a b s should be. Because the multiple scattering between gluons and the pair 
changes the value of Ak c5 while the pair exits the nuclear medium, the effective cr a b s cannot 
be a constant, and it should increase as the cc passes through the medium. Therefore, we 
expect a stronger suppression than that predeicted by the Glauber formula in Eq. ( |3~TD for 
a large nuclear medium, which results into a non-linear dependendence of the suppression 
on the semi-log plot [BH]. The quantitative prediction of the suppression depends on the 



increase of relative momentum per unit length for the cc pairs in nuclear medium. With 
only two parameters: e 2 , the average relative momentum square per unit length in medium 
acquired by the cc pair, and ap which determines the transition probability for a cc pair 



of relative momentum Ak c5 to form a bound-state J/ip meson, it was shown in Ref. |38 
that this picture of J/if> suppression is consistent with all existing data including the NA50 
data on the total J /if) cross section as well as its spectrum (or dependence on transverse 
energy Et of the J/if> events). 

However, in this paper, we are interested in the general features of the aj/^(A,qT) in 
hadron- nucleus collisions at fixed target energies. Because of the fact that the ratio of J/i/j 
total cross sections in hadron-nucleus collisions at fixed target energies can be fitted by the 
Glauber formula in Eq. ([37]), we will not present any detailed model calculation for the ratio 
of J /if) cross sections here ||38|| . Instead, we will approximate the J /if) suppression by using 
the Rjm defined in Eq. ([37]) for our discussions in the rest of this paper. 

Similar to the Drell-Yan process, the nuclear broadening of transverse momentum square 
for J/ip production, A(gf.)j/^, can be calculated in principle within QCD perturbation the- 
ory. In our picture of the J/ip suppresion, the A(q^)j/^ should not be directly tied to the 
J /if) suppression, because the broadening is an effect on the total momentum of the cc pairs 
while the suppression is a result of changing the relative momentum of the pairs. This is 
why the observed nuclear broadening of the transverse momentum square for the J/i/j and 
Drell-Yan production shows a similar A 1 / 3 dependence, while the J/ip and Drell-Yan total 
cross sections have a very different suppression. 

However, due to extra final-state interactions in J /if) production and the difference in 
partonic subprocesses between the Drell-Yan and J/ip production, A(g|>)j/^ will depend on 
the gluon-gluon correlation function in addition to the quark-gluon correlation functions 
mentioned in previous sections |Tj|. But, for a normal nuclear target, both these correlation 
functions should be proportional to the target length or A 1 / 3 . Similar to Eq. (|17]), for the 
normal nuclear targets, we can define 

A(q 2 T ),^ = b,^A l l 3 , (38) 



where can be extracted from existing data. From Ref. [40], bj/^p ~ 0.06 GeV 2 , which 
is larger than bpy- With the averaged transverse momentum square (q^jfip ~ 1-68 GeV 2 
extracted from the fit in Fig. |], we have xj/v> ranges from 0.08 to 0.22 for most relevant 
nuclear targets. Clearly, xj/</> is larger than Xdy, but, it is still a small number for most 
targets. 

Substituting Eqs. ( ]3"7| ) and ( ]3"%| ) into Eq. (]3^|), and taking the small limit, we derive 
the leading contribution to the aj/^(A, qx) for J/if> production in a normal nuclear medium 
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as 



a J/7 p(A,q T ) 



Pr + 



9 J/V 

,2\hN 
T/J/-0 



-1 + 



It 

J2 \ hN 



(39) 



and r is a proportional constant defined as La = r A 1 / 3 . 

l !/, \ rj^ e ] eac }j n g contribution 



we also neglected the very small r(g^)j^ 



where (3 was introduced in Eq. (|37|) 
In obtaining Eq. (|39[ 

to aj/^(A,q T ) given in Eq. ( |3"9"D has the same universal quadratic dependence on q T as 
cxdy(A, qx) in Eq. (p2|) , except a gr-independent shift in magnitude given by — Since 

Lyi is proportional to A 1//3 and A 1 / 3 « hi(A), the shift in magnitude is insensitive to the 
atomic number A. 

Similar to the Drell-Yan case, we obtain from Eq. (BB), 



aj/V>(<?T = 0) 



ln(l + Xj/v) 



\n(A) 



+ 



ln(A) 



In 



a 



hA 



/A' 



a 



hN 

J/V> 



(3r 



(40) 



and it shows that the aj/,/,(gr = 0) is insensitive to the atomic number A for normal nuclear 
targets. Corresponding to ctj/ip(q T ) = 1, the critical value q T for J /if) production is given by 



(It? = (4) 



J/V 



(4) 



1 + 



ln(l + Xj/</>) - In 



a 



hA 

J/V 



/A 



a 



hN 

J/V 



<?t> 



(41) 



which is again insensitive to the atomic number A. 

In Fig. |6], we plot the aj/^,(gr), defined in Eq. fl35|), as a function of g-r for different nuclear 
targets. In plotting Fig. |a, we used (gf.)J% = 1.68 GeV 2 , and A(g|)j / ^ = 0.06A 1 / 3 GeV 2 . 
The solid and dashed lines correspond to the nuclear target C (A=12) and W (v4=184), 
respectively. As expected from Eq. (p9|), Fig. ||a demonstrates that aj/^A, g-r) for J/ip 
production has the same universal functional form as any (A, qx), and it also has scaling on 
A. 

In Fig. |B|a, the thin lines correspond to the large qT region (g-p > 2 GeV), where we 
expect that perturbative calculations of the nuclear dependence should be valid. However, 
such calculations are not available yet. 

Since the suppression of J/ip total cross sections in hadron- nucleus collisions corresponds 
to a g T -independent shift in the value of aj/^(A,q T ), as we explained earlier, the potential 
nuclear dependence for the shape of aj/^(A, q T ) is most sensitive to the nuclear broadening 
A(q T )j/ 1 p. Although QCD calculation of the nuclear broadening of J /if) transverse momen- 
tum square is not available yet, similar to the Drell-Yan case, we expect the broadening 
A(gfi)jM to be proportional to the "soft-hard" parton-parton correlation functions. Differ- 
ent from the "double-hard" parton-parton correlation functions, both position variables for 
soft gluons, such as the efa/f and dy^ in Eq. (|I~4"|), have no exponential oscillation factors. 
Therefore, for a color deconfined nuclear medium, both position variables can be as large 
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as the size of the medium; and in principle, the nuclear broadening for J/ip production can 
be enhanced as much as A 2 / 3 , or A(g^)j/^ rs bA 2 ^ 3 , if we let the color correlation length be 
proportional to A 1 / 3 . Notice, the proportional parameter b given here is not necessary the 
same as the bj/^ extracted from data on a normal nuclear medium. Since we are mainly in- 
terested in the A-dependence of the a (A, gr), we will approximate the b = bj/^ for following 
discussions. 

In a color deconfined medium, the Xi/ip m Eq. ( |3"5] ) is no longer a small number, and its 
J 4 2 / 3 -dependence can not be canceled by the ln(A). Consequently, the aj/^(A, q T ) should 
become very sensitive to the atomic weight A (or the medium size). In Fig. ^|b, we plotted 
the aj/^(A,q T ) with the same parameters as those used to plot Fig. ||a, except A(gf.)jM = 
O.O6A 2 / 3 GeV 2 , which mimics a color deconfined nuclear medium. Notice that in plotting 
Fig. |6|b, we did not include a possible change in A-dependence of the ratio of J/xjj total cross 
sections, since it contributes only to a g^-independent shift in magnitude. From Fig. || it 
is clear that A-dependence of the aj/^iqr) for J/ip production can be a sensitive probe for 
the color deconfinement of the nuclear medium. 

V. DISCUSSIONS AND SUMMARY 

In this section, we discuss the predicting power of the nuclear dependence coefficient 
a(A,qx) for the Drell-Yan and J ftp production derived in the previous two sections, and 
also discuss the uncertainties in the analytical expressions. Finally, we summarize our main 
conclusions. 

The predicting power of the nuclear dependence coefficient a(A, q?) in Eqs. (|2"0| ) and (|3"5]) 
is its universal quadratic dependence on q T , and the fact that all parameters are completely 
fixed by either calculable or independently measurable quantities: (1) the ratio of total cross 
sections, (2) the averaged transverse momentum square in hadron- nucleon collisions, (3) the 
nuclear broadening of the averaged transverse momentum square, and (4) the perturbatively 
calculable r(g|;)'s. The ratio of total cross sections, a hA /a hN , can be either measured in 
experiments or calculated in a theory. The ratio is independent of q?, and therefore, it does 
not affect the shape of a(A, g^). On the other hand, the averaged transverse momentum 
square in hadron- nucleon collisions, (qj-) hN , does not have any nuclear dependence and can 
be independently measured. The nuclear broadening of the averaged transverse momentum 
square, A(gf,), can be factorized into multi-parton correlation functions and perturbatively 
calculable partonic parts. These non-perturbative multi-parton correlation functions are 
universal and can be extracted from the nuclear dependence in the large qT region. Therefore, 
adding the small corrections from the perturbatively calculable r(g^)'s, we can actually 
predict the a(A,qT) for the Drell-Yan and J/ip production. Fig. f|b is an example of our 
prediction and its comparison with E772 data. 

The main uncertainty in our predictions for a(A, q?) comes from the size of the QCD 
resummation for potential large logarithms, such as (a s In 2 (Q 2 / q T )) n . For the Drell-Yan 
production at both collider and fixed target energies, our formula for the o:dy{A, qx) should 
be valid or at least a good approximation, as long as Q 2 is not too large and qr < 2.5 GeV. 
For J /ip production, because Mj/^ is relatively small in comparison with the Q for a typical 
Drell-Yan pair, our formula for aj/^(A, q T ) should be valid for q T < 2 GeV. When q T is 
larger, aj/^(A, qx) should be calculable in perturbative QCD. 
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Another possible uncertainty in our predictions is the asymptotic value of normal parton 
distributions as x — > at low energy, as seen in Eq. (p8|) . The precise asymptotic value of 
the parton distributions as x — > itself is a very interesting quantity, which can provide 



rich information on parton saturation |33fl. For the quantities we discussed in this paper, 
different asymptotic values of the parton distributions result into different normalizations 
of the "double-hard" parton-parton correlation functions. Since the contributions from the 
"double-hard" correlation functions to the o:dy{A, qr) in the large qp region is as important 
as those from the "soft-hard" correlation functions fll|| , different asymptotic values of the 
parton distributions can result into a different crossing point between the low qp spectrum 
and the high qp spectrum, and consequently, a slightly different value of the q^. For example, 
a smaller value of xG(x)\ x ^o corresponds to a slightly larger value of qp, which actually 
describes the data slightly better. Although current data are still not sufficient to provide 
the precise values of parton distributions as x — > 0, such asymptotic values can be measured 
in principle, and corresponding uncertainties can be fixed. 

The four key quantities for determining the a(A,q T ): a hA /a hN , (q T ) hN , A(gf,), and 
r(gy)'s can depend on beam energy, Pbeam, and xf, and also the value of Q in the case of 
the Drell-Yan production. For example, the dependence of a(A, qp) on the beam energy can 
be qualitatively understood as follows. For given values of xp (and Q in the case of the Drell- 
Yan production), a larger Pbeam ( or vS) leads to a larger (q T ) hN because of the larger phase 
space. From Eq. (^) [or Eq. fliPPl, we conclude that when Pbeam increases, the intersection, 
a(qp = 0), with the vertical axis at qp = in Fig. |4] [or Fig. ^|a] moves closer to one. On the 
other hand, from Eq. fl24|) [or Eq. fl4lDl, the increase of the beam energy leads to a larger 
value of the critical point qp at which the a(A, q^) = 1. Therefore, the curve of a(A,qp) 
becomes flatter when Pb ea m increases, which is consistent with the data from NA10 ||41|| . 
Similar qualitative conclusions can be derived for the ^-dependence (or Q-dependence in 
the case of the Drell-Yan production) of a (A, qp). In any case, one should take into account 
the possible xp dependence of all four key quantities that determine a(A, q T ). 

Recent data from Fermilab experiment E866 shows that for different regions of xp, 
the aj/^(A,qp) for J/ip production have similar shapes in g^-dependence, but, different 
magnitudes ||. It was found that for three xp regions (small, intermediate, and large), the 
magnitude of aj/^(A, qp) decreases as xp increases. This phenomenon is consistent with our 
prediction of aj/^(A, qp) in Eq. fl35|) [or Eq. (|39|)1. According to Eq. (^), the suppression in 
J/ip total cross sections represents a gr-hidependent shift in the magnitude of aj/^(A, qp). 
The larger the suppression, the smaller aj/^(A,qp). Because of limited xp range, the ratio 
of total cross sections, Rf/^ = (V^) a j / v/ cr J/v> ^ n ^l- (HD> becomes xp dependent. It is 
an experimental fact that the larger xp is, the more suppression for J/?/> production (or 
smaller Rf/^) [0. Therefore, from Eq. fl35[), the larger xf, the smaller aj/^(A, qp), which 
is consistent with experimental data ||. Although we do not have the four key quantities 
for different xp regions to make absolute predictions, we can still test the universality (or 
consistency) of the aj/^(A,q T ) by using the quadratic form in q T shown in Eq. (^) to fit 
the data. In Fig. 0, we plot E866 data on aj/^(A,q T ) in three separate xp regions: small 
(SXF), intermediate (IXF), and large (LXF). We also plot in Fig. [7| our fits to these data 
with a universal quadratic form in qp, as shown in Eq. fl35|). It is clear from Fig. [7] that 
our universal function for aj/^(qp) is consistent with all data in the small qp region, which 
covers qp < qp ~ Mj/^/2. Also, as expected, our formula for aj/^(A,q T ) deviates from the 
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data when q? > Mj/^,/2. Like the Drell-Yan case, shown in Fig. |], it will be very interesting 
to calculate the nuclear dependence of the J/ip transverse momentum spectrum in the large 
qx region to test the QCD dynamics. 

In summary, we derived an analytic formula for the nuclear dependence coefficient 
a(A, q?) for both the Drell-Yan and J/i/j production in the small qx region. The formula has 
a universal quadratic dependence on qx, and all parameters are completely determined by 
four key quantities: a hA /cr hN , (q T ) hN , A(q T ), and r(g^)'s. These quantities can be either 
calculated in QCD perturbation theory or independently measured in other experiments. 
We explicitly demonstrated that our a(A,qx) is consistent with existing data. 

Furthermore, we showed that our a (A, q?) is extremely insensitive to the atomic weight 
A for a normal nuclear target. However, for a color deconfmed nuclear medium, the a(A, qx) 
becomes strongly dependent on A (or the medium size). Therefore, the A-dependence of 
a(A, qr) can be a potential observable for detecting a color deconfmed system. 
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FIG. 1. Typical Drell-Yan transverse momentum spectrum as a function of qx in three regions. 
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FIG. 2. Data on the Drell-Yan transverse momentum spectrum: (a) E772 |^] and (b) CFS 
[27], in comparison with a Gaussian-like fit in the small qx region plus a calculated perturbative 
tail in the large qx region. 
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Fig. 3 

FIG. 3. Sample Feynman diagram for the double scattering contributions to the nuclear de- 
pendence of the Drell-Yan gj-spectrum at large qx- 
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Fig.4 (b) 

FIG. 4. (a) a{A,qx) for the Drell-Yan production as a function of qp. Low qT region is given 



by Eq. (|20|), and high qx region is given by Ref. [11| with the maximum and minimum values of 
the C discussed in the text, (b) R(A,qx) for the Drell-Yan production as a function of qx- Data 
from E772 P,|34]|, and the theory curves are given by R(A,qx) = A a<yA,qT ^ 1 with a(A, qx) given in 
(a). 
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Fig. 5 

FIG. 5. Fermilab E771 data on transverse momentum spectrum [35| in comparison with a 
fit of Gaussian-like distribution. 
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FIG. 6. a(A,qx) for J/tp production, denned in Eq. (|3q), as a function of qx for normal color 
confined nuclear targets (a), and for possible color deconfined nuclear medium (b). 
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Fig. 7 

FIG. 7. Comparison of fits, using our universal quadratic form for a(A, qx) in Eq. (^), with 
J/^ data from Fermilab E866 collaboration || in three regions of different xf. 
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